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Abstract. Conformational properties of polymer melts confined between two hard 
structureless walls are investigated by Monte Carlo simulation of the bond-fluctuation 
model. Parallel and perpendicular components of chain extension, bond-bond 
correlation function and structure factor are computed and compared with recent 
theoretical approaches attempting to go beyond Flory's and Silberberg's hypotheses. 
We demonstrate that for ultrathin films where the thickness, H, is smaller than the 
excluded volume screening length (blob size), ^, the chain size parallel to the walls 
diverges logarithmically, /2N « &^ -|- clog(A^) with c ^ 1/H. The corresponding 
bond-bond correlation function decreases like a power law, C{s) = dj with s being 
the curvilinear distance between bonds and uj = \. Upon increasing the film thickness, 
if, we find - in contrast to Flory's hypothesis - the bulk exponent w = 3/2 and, 
more importantly, a decreasing d{H) that gives direct evidence for an enhanced self- 
interaction of chain segments reflected at the walls. Systematic deviations from the 
Kratky plateau as a function of H are found for the single chain form factor parallel 
to the walls in agreement with the non-monotonous behaviour predicted by theory. 
This structure in the Kratky plateau might give rise to an erroneous estimation of the 
chain extension from scattering experiments. For large H the deviations are linear with 
the wave vector, q, but are very weak. In contrast, for ultrathin films, H < ^, very 
strong corrections are found (albeit logarithmic in q) suggesting a possible experimental 
verification of our results. 
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1. Introduction 

The excluded volume interaction or self-repulsion of the segments along an isolated chain 
leads to self-avoiding walk statistics of long polymers in dilute solution [TJ |2l El IH EI ■ 
The chain extension as measured by the mean squared end-to-end distance, R^, increases 
like a power law R ~ N'^saw -^i^h chain length A^. The exponent adopts the non-trivial 
value I'sAw = 0.588 • • •. The chains swell in order to reduce intramolecular interactions 

[HE!. 

In a semi-dilute solution or a melt, however, these excluded volume interactions 
are screened on length scales larger than the size of the "blob", ^ pi 121 El IE]- For 
distances smaller than ^ the chain statistics obeys self-avoiding statistics while for 
larger distances the chains adopt Gaussian conformations which are characterised by the 
exponent u = 1/2. On large distances the chains can not reduce the segmental repulsion 
by increasing their spatial extension, the majority of excluded volume interactions stems 
from intermolecular contacts p]. 

The concept of screening of excluded volume interactions in a dense polymer melt 
- the Flory hypothesis - plays a pivotal role in modern theory of polymer melts and 
concentrated solutions. It has been corroborated by Edwards jSl IH| and is born out 
of renormalisation group calculations |6|. The description of chain conformations as 
Gaussian random walks lays at the heart of many analytical or numerical approaches 
like, for instance, the self-consistent field theory ^] which has been widely used 
to describe spatially inhomogeneous polymer systems like surfaces of polymer melts, 
interfaces in polymer blends or self-assembly of copolymer systems. 

One consequence of screening of excluded volume interactions on length scales larger 
than ^ is the behaviour in a thin film (see Fig. ^^j: If the chains were describable by 
random walks without long-ranged interactions the properties parallel and perpendicular 
to the film surfaces would decouple. According to Silberberg's hypothesis [l2j the chain 
conformations could be conceived as random walks reflected at the film surfaces, i.e., 
parallel chain extensions would remain unperturbed. 

As we decrease the film thickness, H, the chain folds back into the volume occupied 
by their coil and other chains are expelled from that region. Silberberg's hypothesis 
has to break down if the self-density is comparable to the segmental number density 
of the bulk melt, p, i.e., pR^H ~ N Olllj. For H <^ the chains adopt quasi-two- 
dimensional configurations j2] and the number of intermolecular contacts are strongly 
reduced. This is illustrated in Fig. mb). 

While these phenomenological considerations yield a description of the leading 
scaling behaviour of large scale chain conformations, there are subtle corrections that 
may lead to significant deviations: Renormalisation group calculations provide great 
insight into the way excluded volume interactions are screened upon increasing the 
density and crossing over from dilute to semi-dilute solutions [2111 Ej- These calculations 
predict the blob size, ^, and the chain extension, R, as a function of the strength 
of the excluded volume interaction and the segmental density, p. They also provide 



Thin polymer films 



3 



information about the crossover from self-avoiding walk behaviour to Gaussian statistics 
upon increasing the distance along the chain in semi-dilute solutions. Surprisingly, 
renormalisation group calculations by Schafer [Bj and Monte Carlo (MC) simulations 

HSI of the bond-fluctuation model revealed deviations from the plateau in the 
Kratky plot q'^S{q) vs. q which were traced back to a Goldstone-type singularity. These 
deviations from the plateau in the Kratky plot indicated deviations from the Gaussian 
behaviour on length scales /anjferthan ^. Upon confining the chains into a thin film, one 
observes an increase of the deviations from the Gaussian behaviour fL3\, IT^ . 

Recently, the statistics of quasi-two-dimensional polymer melts (called "self-avoiding 
trails") has been studied by Semenov and Johner ^l]. This model corresponds to 
polymers confined to ultrathin films, H <^ ^, where, however, chains and chain segments 
are still allowed to overlap as opposed to the pure two-dimensional limit {H 0) which 
has been considered thoroughly in the past resulting in beautiful analytical predictions 
[in] and more recent numerically tests [T71 UHl UHl 120] • Rather than calculating the 
chain extension as a function of density and strength of the excluded volume interaction 
Semenov and Johner consider the residual excluded volume interaction which is a weak 
perturbation to Gaussianity in D = 3, but is marginal, i.e., has to be renormalised, in 
ultrathin layers. This allows for quantitative predictions in the melt regime both for 
films and in the bulk. In contrast to Silberberg's argument, the chain extension for these 
D = 2 + e systems is predicted to exhibit logarithmic corrections, 

RliN)/2N = b{H)^ + ciH)\ogN (1) 

for the parallel components {a = x, y) of the end-to-end distance with respect to the 
walls. The first coefficient, b, corresponding to the statistical segment length is predicted 
to depend only very weakly on H [TT]. The second coefficient c{H) should be positive and 
inversely proportional to the number of particles per unit surface and, hence, inversely 
proportional to the ffim thickness, H. Likewise, one finds in the bulk [D = 3): 

Rl{N)/2N = h\p)-c{p)/^ (2) 

for any component a = x,y,z of the chain end-to-end vector [21]. Note that the 
second coefficient decreases now inversely with the monomer density. In both cases, 
the coefficient c of the leading correction term does not depend on the excluded 
volume parameter, v, expressing the fact that the corrections are due to the large-scale 
incompressibility of the polymer melt [Tl| l2T] . 

In our manuscript we numerically investigate the chain statistics in thin polymer 
ffims paying particular attention to residual excluded volume effects and their 
dependence on the ffim thickness. As sketched in Fig. ^ our computational study 
covers the regime from thick films, H > R, to ultrathin films, H < ^, and allows 
us to investigate the crossover from the bulk behaviour {D = 3) over ultrathin films 
(D = 2 + e) up to purely two-dimensional systems {D = 2). In the next section we 
describe the model and simulation technique. Then we investigate the chain extensions 
parallel and perpendicular to the ffim, the bond-bond correlation function and the single 
chain structure factor. The final section presents our conclusion. 
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(a) D=3 : R, H » ^ (b) D=2+e : 




Figure 1. Schematic illustration of a polymer chain with typical size R confined 
between two parallel hard and structureless walls at a distance H. The z-axis is 
perpendicular to the lower wall. Periodic boundary conditions are used in x and y 
direction, (a) When ^ ^ R <^ H, being the excluded volume blob size, coils far 
from the walls are unperturbed, i.e., Ra — i?buik for all components a — x,y,z. 
According to Flory and Edwards, they can be represented, to first approximation, as 
Gaussian chains of blobs of size ^. For chains (or chain sections) close to the walls it 
is generally assumed that they may be reflected at the walls but remain Gaussian, (b) 
For 7J ~ ^ <C i?buik, the coils become effectively two-dimensional, but are still allowed 
to overlap and to cross in the two-dimensional projection. These are called D = 2 + e 
systems in contrast to a purely D — 2 films where intersections are forbidden and the 
coils become compact. It is illustrated that in the D = 2 + e case, the correlation 
length parallel to the film surfaces differs from H. 



2. Computational model and some technicalities 

Coarse-grained models are very efficient to investigate the universal properties of dense 
polymer systems. In these kind of models, one integrates over the microscopic degrees 
of freedom, which do not affect the universal physical laws but only their prefactors 
and local properties. In this work we have used the "bond fluctuation model" (BFM) 
- a lattice Monte Carlo scheme where a monomer occupies 8 lattice sites (i.e., the 
volume fraction = 8p, p being the number density) and the bonds, 1, between adjacent 
monomers can vary in length and direction, subject only to excluded volume constraints 
and entanglement restrictions j221l23- The BFM has been extensively studied in the 

past pnaiiHiiniiiHiETiEniEn!. 

Chain configurations are updated by two types of canonical moves: random 
monomer hoppings and slithering-snake (or reptation) moves [21] • In the first case 
we consider random displacements of an effective monomer by one lattice site in one of 
the 6 possible directions in the lattice. These moves are efficient in relaxing the local 
properties of the chain, such as the bond angle or the bond length. In slithering-snake 
moves, a segment of the chain is removed from one end and added to the opposite end 
of the chain in a random direction. The latter moves relax the chain conformation a 
factor N faster than the random monomer hoppings. (See ref. jSnj for details on the 
exponential increase of the relaxation time for very long chains due to correlations in 
the slithering-snake motion.) 
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As illustrated in Fig. d we use a L x L x H box with hard walls at z = and 
z = H + 1, i.e., the thinnest film with H = 2 corresponds to a pure D = 2 system where 
chains are not allowed to overlap. For H = 4 this becomes possible and, therefore, this 
is an ideal system to test the predictions of ref . [TT] . We have systematically varied the 
film thickness between H = 2 up to H = 84 while keeping the parallel dimensions to 
the walls fixed at L = 256 (with the exception L{H = 2) = 512). Periodic boundary 
conditions have been used in x and y directions. Note that L is always much larger than 
Rx, the typical chain size component parallel to the walls. 

The chain length ranged between N = 16 up to 256. We have concentrated in this 
study on dense polymer melts and all the data reported here are for a number density 
p = 0.5/8. For this density the (bulk) excluded volume length is about ^ «i 7 

We used the configurational bias method to create the initial configurations for 
chain lengths up to = 64. (The longest chain length possible numerically with this 
method decreases with H.) To this end a chain is grown monomer by monomer into 
the system taking into account all possible bond vectors for the next monomer. From 
those positions of the next monomer one that does not overlap with any other monomer 
in the system, is chosen at random. This choice introduces a bias which is removed 
by the Metropolis acceptance criterion once the chain is fully grown. For larger A^, 
we started with relatively compact two-dimensional coils {R ~ N^^'^) oriented parallel 
to the surfaces. In all cases we equilibrate our systems using a mixture of local and 
slithering snake moves which has been found to be the most efficient Some of 

our configurations have already been used in some recent related investigations on thin 
polymer films ^3 • 

Various static and dynamical properties parallel and perpendicular to the walls 
have been measured on the fiy and for a more detailed analysis configurations were 
periodically stored. Since the x and y directions parallel to the walls are equivalent these 
properties are averaged together. If not specified otherwise, properties are averaged over 
all chains irrespective of their distance to the walls. For later reference, we note here the 
bulk values for the end-to-end distance R ~ 49.8 and the radius of gyration Rg ^ 20.3 
for chains of length A^ = 256. The three components of the mean-squared bond vector 
= (^2)^/2 _ ^_g2 of the bulk remain unchanged for all H > 7. Obviously, the 
perpendicular component must ultimately vanish and the parallel component increase. 
For if = 4 we find, for example, Ix = ly = 1-63 and 1^ = 1.2. 

3. Results 

We discuss now in turn three intra-chain properties: the global chain size Ra{N, H), the 
bond-bond correlation function Ca{s) as a function of curvilinear distance s along the 
molecule's backbone and the form factor Sa{q)- 
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3.1. End-to-end distance and radius of gyration 

We start the discussion by presenting the chain size as measured by the three components 
i?^ = {{ra,N — f^a,iY) of the end-to-end distance R and the corresponding components 
of the radius of gyration Rg. Here, r^^n denotes the a-component of the position of the 
nth monomer of the chain. 

In Fig 121 we present the end-to-end distance and the radius of gyration, Rg, as 
a function of the film thickness H. The data are normahsed by the corresponding 
bulk value, i.e., all the ratios must become unity for large H. Only data for chain 
length N = 256 has been included. As expected, the thinner the film, the larger 
is the component parallel to the surface and the smaller is the perpendicular one. 
The perpendicular component increases continuously with H. The effect is linear 
for H <^ -Rbuik, as one expects even for perfectly Gaussian chains. Note that the 
perpendicular component of Rg increases more rapidly as the component of the end- 
to-end vector for small H. In agreement with ^l] the (reduced) parallel components 
remain constant for H ^ H* ^ ,^ ^ 7, but increase rapidly for thinner films. (We remind 
that for if = 2 no chain overlap is possible in our model and the chains are rigorously 
two-dimensional.) This demonstrates that perpendicular and parallel components of 
the polymer couple in contrast to Silberberg's hypothesis. Note that the observed chain 
length dependence of the crossover slit width H* is relatively weak for the chain lengths 
we have probed (cf. Fig. |2l inset). This is in line with the predicted logarithmic increase 
H* ^ ,^ log(A^/5'), ~ denoting the numbers of monomers contained in the excluded 
volume blob, which follows readily from the condition pR^H* ~ utilising Eqs. (0) 
and (HdH. 

In the next figure we present the parallel component of the end-to-end vector, R^, 
as a function of for different H. Also included is the data for bulk systems {D = 3) 
without walls and periodic boundary conditions in all directions. In agreement with 
Eq. we find corrections to ideality even in this latter case. For details the reader 
is referred to the recent study |2I]. As expected, the reduced data of both D = 3 
and D = 2 limits {H = 2) become chain length independent for long enough chains. 
(Note that the chain size in D = 2 is R^ ~ A^^''^ without logarithmic correction [TE].) 
More interestingly, we find that all intermediate data diverge logarithmically for long 
enough chains (thin lines). Extrapolating our data, one even expects that for A^ ^ 2000 
and H = 4 the chains become actually larger than the pure D = 2 ones. In their 
work Semenov and Johner ^l] have predicted for the slope in log-linear coordinates 
[cf. Eq. ©] 

l/c = Attco = AnpH (3) 

where cq is the number of particles per unit surface. Hence, the coefficient decreases 
inversely with H. In contrast, the first coefficient, b, of Eq. ((H) is predicted to depend 
only very weakly on H Using Eq. Q and adjusting b{H) we obtain rather good 
agreement with the numerical data. This is the central result of this work. Note that 
b{H) is essentially constant as can be seen from the convergence of the lines for small A^. 
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Figure 2. Parallel and perpendicular components of end-to-end distance R (closed 
symbols) and of radius of gyration Rg (open symbols) for TV = 256 as a function of the 
film thickness H. All the quantities are normalised with respect to the corresponding 
bulk value. The data are averaged over all chains in the film irrespective of the distance 
of their centre of mass from the walls. The (reduced) parallel components remain 
constant for distances larger than the characteristic width H* which for the chain 
length used is of the order of ^. The bulk value of ^ for the presented volume fraction 
is indicated. Note that ioi H = 2 no chain overlap is possible in our model and the 
chains are rigorously two-dimensional. D = 2 + e behaviour is expected for 2 < < ^. 
The perpendicular component increases continuously with H . The effect is linear for 
H <C -Rbuik, as one expects. The inset presents the chain length dependence of the 
parallel chain dimensions as a function of H/H* for different chain lengths, TV, as 
indicated in the key. 



We even get acceptable fits for film widths three times larger than ^ which is consistent 
with the logarithmic increase of H* with chain length mentioned above. 

We have also investigated the components of the mean-squared distance along the 
chain -Ra(s) = ((ra,n — '"a,n-fs)^)„ obtained by averaging over all chains and all pairs of 
monomers separated by a curvilinear distance s = \m — n\ along a chain. Here, (■ ■ ■)„ 
denotes the average over all couples of bonds having curvilinear distance s and over all 
polymers in the system. The predictions of ref. ^HElj, Eq. ((1} and Eq. (j2)), generalise 
readily from s = iV to arbitrary curvilinear distance. This is well confirmed by our data 
which looks quite similar to the one given in Fig. |H1 Instead of discussing these data we 
will rather proceed by discussing the numerically more challenging curvature of Ra{s), 
i.e., by presenting our results on the bond-bond correlation function, which allows a 
much more accurate test of the predictions. 
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Figure 3. Parallel components of the end-to-end distance for different H a.s a function 
of chain length N. We plot R^/2N in log- linear coordinates to demonstrate the 
expected logarithmic divergence for Z3 = 2 + e systems The thin lines compare 
with the prediction, Eqs. Ij 113(1 . where we directly verify c{H) = l/AirpH and fit for 
the coefficient b{H) which is found to depend very weakly on the film width becoming 
ultimately H independent (left side of figure). The pure D — 2 data {H ~ 2, stars) 
become chain length independent for chain lengths > 64 (dashed line). Note that 
even the bulk data is monotonously increasing in agreement with Eq. 



3.2. Bond-bond- correlation function 

The connection between mean-squared distance between two segments n and m along 
the chain and the bond-bond correlation function is due to the exact formula 

1 (9^ 

{la,nla,m) = — " Q^Q^ {{^a,n — fa,m) ) (4) 

valid for a = x,y, z. Note that the bond vector is defined as In = r„+i — r„ ^ drn/dn., 
the latter identity becoming true in the continuous limit. More generally, Eq. Q may 
be written in terms of the difference operator A„/(n) = f{n -|- 1) — f{n). In that sense, 
it remains valid for arbitrary curves, r„, even for non-differentiable Brownian walks. 

Averaging again over all pairs (ra, m = n + s) we have computed the correlation 
functions Ca{s) = {la,nla,n+s)n /^a- § Neglecting chain end effects and using Eq. (jH) one 
obtains the compact relation 2Cq,(s)/^ = (9f-R^(s) between both curvilinear properties. 

§ We have compared this definition to {ea.ne-a,n+s)n ^r the normalised bond vector e„ without finding 
any measurable difference. This might be due to the weak correlation between bond length and bond 
angle fluctuations. 
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Figure 4. Perpendicular bond-bond correlation function, Cz(s), vs. curvilinear 
distance, s, between monomers for N — 256 and different film thicknesses, H. The 
inset shows the behaviour of the same correlation function for polymers with centre of 
mass at distances z = 4 or 30 from the wall at z = in a film with H = 62. 



From Eqs. (P) and (0), generalised to arbitrary s, we immediately see that the 
bond-bond correlation function must become 

C„(s) = ds-^ (5) 

with a pref actor ~ 1/p and an exponent u = 3/2 for the bulk and d = c/l"^, u = 1 for 
ultrathin films. Noticeably, the long-distance behaviour of Ca{s) is not characterised by 
a length scale (persistence length) but the residual excluded volume interaction change 
the dominant scaling behaviour to a scale-free power law. This is in marked contrast to 
Flory hypothesis where long-ranged correlations between bonds are neglected and Cq,(s) 
is assumed to decay exponentially (or even faster) |2I]. The technical advantage of the 
bond-bond correlation function compared with the curvilinear distance is that it does 
not depend on the trivial, but very large Gaussian contribution which drops out after 
differentiating. Hence, it allows us to focus directly on the corrections to the Gaussian 
behaviour on large length scales. 

In Fig. m the perpendicular bond-bond correlation function, Cz{s), is plotted as 
a function of the curvilinear distance, s, for chain length = 256 and different film 
thicknesses, H. For very thin films negative correlations at short distances can be 
observed. For H = 4 and H = 5, on average, the angle between two consecutive 
bonds, in the direction perpendicular to the walls, is larger than 90° (i.e., the two 
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vectors have opposite direction). For H = 7, the correlations reach the minimum for 
s = 2. This means that chains fold back every two monomers. These results are a 
trivial consequence of confinement. For very thin films, due to strong fluctuations of the 
density as a function of the distance from the walls, the system essentially consists in the 
superposition of mono-layers with distance equal to the lattice spacing. For H = 4 and 
H = 5we have two layers, then two consecutive segments along the chain are prevalently 
located in different planes and the chain is reflected at the surface. For H = 7, we have 
one more layer and then up two three linked monomers can be placed at different height. 
Upon increasing the film thickness, H, above the blob size, ^, Cz{s) gradually tends to 
the curve obtained in the bulk. 

The same effect can also be observed for thicker films taking into account only 
chains with centre of mass having a distance, z, from the walls smaller than ^. This is a 
consequence of excluded volume effects in the direction perpendicular to the wall which 
give rise to a strong anisotropic deformation of the coil. In the inset of Fig. |3] we show 
Cz{s) calculated for polymers at different distance, z, of their center of mass from the 
surface. In this plot we compare the curves obtained for z = 4 and z = 30 in a film with 
if = 62. While in the first case the interactions with the wall must be very strong the 
walls should have a negligible effect on the latter where the chains are exactly located 
in the middle of the film. Indeed, only in the first case evident negative correlations can 
be seen and we recover bulk behaviour in the second. The insufficient statistics does 
not allow a precise location of the minimum for z = 4. This value roughly corresponds 
to s = 5 which is close to the distance of the center of mass from the wall. 

A quantitative analysis of the decay of C^i^s) is given in Fig. El where we plot the 
function on a double logarithmic scale. The curves can be fitted rather well by power 
laws, in agreement with Eq. Q, and exponential behaviour can definitely be ruled 
out for long enough chains. Obviously, this is a clear-cut contradiction to Flory's and 
Silberberg's descriptions. For large thicknesses, if > 7 ~ ^, we observe the expected 
bulk exponent u = 3/2. Interestingly, the prefactor, d, increases as we decrease the 
film thickness, H, although the exponent remains constant. This shift can be related 
to the dependence on the fraction of polymers oc 1/H which interact with the surface. 
Chain segments close to the wall have increased bond-bond correlations {lx,nlx,n+s) due 
to the enhanced self-interaction of the segment which is, in turn, due to the reflection 
of the segment at the wall. To first order, the reflections may still be described by 
Gaussian statistics. As H increases and the relative population of chain segments of 
size s decreases the amplitude must decrease. This interpretation has been directly 
confirmed by computing, as in Fig. HI Cx{s) for various distances, z, of the chain center 
of mass from the walls (not shown). For ultrathin D = 2 + e films, the bulk exponent 
must break down. Indeed, our data are still compatible with a power-law but with the 
uj = 1 exponent predicted for this regime. 

Since finite-chain size effects could be crucial here larger chains are needed to 
rigorously establish the exponent over more than a decade in s. In the bulk the 
correlation function decays more rapidly for short chains and its explicit dependence 




s = |m-n| 



Figure 5. (a) Log-log plot of Cx{s) vs. s for different film thicknesses and for N — 256. 
The lines indicate the power-law behaviour, Cx{s) = d{H)/s'^, predicted for different 
regimes. The slope cj = 3/2 (bottom line) is expected for bulk systems and larger 
slits. While this exponent remains unchanged for H ^ ^ the amplitude d{H) does. 
As indicated by the broken line, we demonstrate a systematic increase with decreasing 
H. Finally, for I? = 2 + e films our data confirms the exponent uj ~ 1 (upper bold 
line) implicit to the work of Semenov and Johner jjl]. (b) Double logarithmic plot of 
Ccnd,a:(s) in thc bulk and a thin film, H — 4, corresponding to = 2 -|- e. The bulk 
correlation function averaged over all pairs is shown as reference. The predicted power 
laws are indicated by straight lines. 

on has been obtained [25- This is, in part, due to the behaviour close to the 
chain ends. In the analysis we have averaged the correlations over all pairs with the 
same distance, s, along the chain's contour but neglected chain end effects. This is 
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permissible for large A^.|| Interestingly, the correlation function measured from a chain 
end, Cend,a{s) = {la,ola,s)n /^a decays with a different, stronger exponent, c^e = -D/2 + 1. 
The simulation results are presented in panel (b) of Fig. El For bulk melts, very long 
chains are available [21] and confirm the value of the exponent, = 5/2 in D = 3. The 
data for thin films are also compatible with the theoretical prediction. 

3.3. Form factor S{q) 

The single chain form factor 

... 1 A / ..^ . .A 



of a chain is an important, experimentally relevant quantity. The components of the 
scattered radiation vector are denoted by qa- We are only interested here at the internal 
correlations parallel to the wall and, hence, consider only wave vectors parallel to the 
surface {qz = 0) and omit the index in the following. 

In Fig. iniwe present S{q) for different film thicknesses and chain length N = 256. 
The analog quantity, calculated for the bulk system, is also depicted as a reference 
curve. For small q the form factor just counts the number of scattering units, for large 
q it probes the structure on the scale of the monomers, the so-called Bragg peak. Not 
surprisingly, in both limits S{q) does not depend on the presence of walls. For small q the 
data should be described by the rather general expansion S{q) = N[l — {Rg^qY + ■ ■ ■] 
due to Guinier Here, Rgx{H) stands for the measured parallel component of the 
radius of gyration. Since Rg^ increases for small H (Fig. 12)) the structure factor should 
decrease in agreement with what can be seen in the figure. One can check explicitly for 
all H that Guinier's formula applies. Hence, the deviation from the bulk limit seen in 
this regime is fully described in terms of the increased chain size. 

The Gaussian self-similar structure for intermediate length scales is the essence of 
Flory's hypothesis and results in a power-law behaviour of the form, S{q) ~ l/q"^, which 
is indicated by the broken line. Broadly, this is confirmed by the thick films {H > 
Differences become more apparent when we plot the ratio of the structure factor of a 
Gaussian chain, S'nebye, with the same radius of gyration and the measured single chain 
structure factor. This ratio is presented in the inset of Fig. IHl The clearly observable 
deviations are already present in the bulk ^Hj and increase upon reducing the film 
thickness |lH]. Significant differences with respect to Gaussianity are observed for the 
thin films. 

The deviations from the Gaussian behaviour as well as the effect of the confinement 
are magnified in the Kratky plot shown in Fig. [7| where we have plotted {qRgxyS{q)/N 
vs. qRgx for = 256. This plot illustrates a popular method for extracting the chain 
dimensions from scattering data. Using the measured value of the radius of gyration we 
obtain perfect scaling for the Guinier regime, qRgx < 1- Not surprisingly, the scaling fails 

II The dependence of correlations on the position along the chain has recently been studied for isolated, 
self-repelling polymer chains by Schafer and Eisner PB] 
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Figure 6. Single chain form factor for wave vectors q parallel to the wall surfaces 
for different H. Only data for N ~ 256 and p — 0.5/8 are considered. The bulk 
result {H oo) is given by the bold line, Flory's hypothesis indicated by the broken 
line. The Bragg limit {q > 1) is not affected by the presence of the walls. However, 
this can be clearly seen for larger length scales where the scattering amplitude goes 
systematically down with decreasing film width. The effect is the most prominent for 
D — 2 system (stars), but even there the Bragg peak is not altered. The inset presents 
the ratio of the structure factor, ^Debye, of a Gaussian chain with the measured radius 
of gyration and the structure factor, S{q), observed in the simulations as a function of 



for very large q (Bragg peak) where the data diverges. More interestingly, the scaling 
fails even for intermediate wave vectors where it should hold for Gaussian statistics [T5] . 
Instead we find a pronounced non-monotonous behaviour in stark contrast to traditional 
hypotheses. As already seen in the previous figure, the deviations get more pronounced 
with increased confinement ^^l- Similar deviations are also observed in semi-dilute bulk 
solutions of long chains, N = 2048 [15^ . 

These effects have been analytically predicted in semi-dilute solutions and melts 
[TT]. For qRgx ^ 1 and for long chains ^ 



has been suggested for thin films. Only the first of the two coefficients e{H) and f{H) 
should depend sufficiently strongly on the slit width to be measurable: e{H) ~ c{H) ~ 
1/H [TT]. SDehyeil) deuotes the structure factor of a Gaussian chain with identical 

% For exponentially long chains yet another regime is expected |11 | which is outside the reach of 
simulation. 




l = e{H) log(g/(iJ)) 



(7) 
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Figure 7. Kratky plot {qRgx)'^ S {q) / N vs. qRgx for wave vectors parallel to the walls 
using the same symbols as in the previous plot. We have rescaled the data with the 
chain length N and, Rgx, the measured component of the radius of gyration parallel 
to the surface. We see that the density fluctuations get systematically suppressed with 
decreasing film width H . The data for H — A has been compared with the logarithmic 
correction, Eq. proposed in ref. (llj and the bulk data (bold line) with the 
corresponding linear relationship (dashed line). The dashed-dotted line indicates the 
Debye functions SochyciqRgx) / N . The inset presents the deviation, S'Dcbyc(9) / S{q) — 1, 
of the measured form factor in thin films with 4 < H < 12 from the Debye function 
for the same radius of gyration. As expected, for 1 > q > 1/Rg the deviations grow 
logarithmically in q and the amplitude is proportional to 1/H as can be seen from 
the scaling collapse. The deviation decreases again for large q due to additional (non- 
universal) corrections in the Bragg limit which arc not taken into account by the 
theory. 



radius of gyration. In the inset of Fig. [7| we explicitly verify the film thickness and 
wave vector dependence for 4 < H < 12. It also confirms that e{H) log{f{H)) is rather 
small. For D = 3 melts the deviation from the Kratky plot has been estimated more 
recently by Beckrich |2Z] using the perturbation approach of ref. ^T]. This yields the 
linear expression ^^^g^^^-^^'^^ — 1 = e(p)g which we have fitted to our bulk data (dashed 
line). A similar expression was derived for the limit of semi-dilute solutions (p ^ 1 but 
pR^/N > 1) by Schafer S'Debye(g)/5'(g) - 1 = e'[(gOT"'''^ + ••• where e' is a constant 
and e = 4 — D [6 . Thus we expect corrections to the Kratky plot to be present over 
the entire concentration regime where the chains overlap. Since not all pref actors are 
explicitly available as a function of the model parameters of the BFM the test is less 
rigorous than the one in Fig. El While our simulation data are compatible with the 
theoretical predictions longer chain lengths would clearly be desirable. 
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4. Conclusions 

Using Monte Carlo simulations of long polymers confined between hard structureless 
walls we investigate deviations from the classical hypotheses of Flory P and Silberberg 
[T^ . In agreement with recent theoretical work JT] on ultrathin films with if < ^ we 
find that parallel and perpendicular components couple and the chain size parallel to the 
surfaces increase strongly. We quantitatively demonstrate the logarithmic divergence, 
Eq. for the chain size parallel to the walls as a function of chain length, N. The 
prefactor c ^ 1/H predicted for the logarithmic deviation allowed a rather accurate 
fit of our data. As emphasised in refs. [HI |2l], these effects are due to the chain 
connectivity and the large-scale osmotic incompressibility of the solution and express 
universal physics independent of local properties. 

Perhaps the most striking effect of this study is the power law asymptote for the 
bond-bond correlation function which measures the curvatures of the curvilinear distance 
along the chains. It allows a direct numerical test of the non-Gaussian corrections both 
in the bulk and in the ultrathin films and demonstrates the presence of long-ranged 
correlations neglected by the classical hypotheses. Note that the decay exponent for 
both bulk and thin film geometry can be expressed in the form u = uD, D being 
the effective dimensionality of the system and u = 1/2 the characteristic exponent of 
a random walk. Interestingly, this is similar to the result obtained for the velocity- 
correlation function of dense liquids in two and three dimensions fIE\ . 

An important consequence of our work arises for an experimentally relevant 
quantity, the static structure factor S{q). In fact, simulation and theory show distinct 
non-monotonous behaviour of q^S{q) vs. q (Kratky plot) due to the non-Gaussian 
corrections which even get enhanced with decreasing film width. This suggests a possible 
route for experimental verification and is cause of serious concerns with regard to the 
standard operational definition and measure of the persistence length from the "assumed" 
Kratky plateau. 

Finally, we point out that the physical mechanism which has been sketched above 
is rather general and should not be altered by details such as a finite persistence length 
— at least not as long as nematic ordering remains negligible. This is in fact confirmed 
by preliminary and on-going simulations. 
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